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ABSTRACT 



We present an analysis of the spin and orbital properties of the newly discovered accreting pulsar IGR J 17480-2446, located in the 
globular cluster Terzan 5. Considering the pulses detected by the Rossi X-ray Timing Explorer at a period of 90.539645(2) ms, we 
derive a solution for the 21.27454(8) hr binary system. The binary mass function is estimated to be 0.021275(5) M Q , indicating a 
companion star with a mass larger than 0.4 M Q . The X-ray pulsar spins up while accreting at a rate of between 1.2 and 1.7xl0" 12 Hz 
s -1 , in agreement with the accretion of disc matter angular momentum given the observed luminosity. 

We also report the detection of pulsations at the spin period of the source during a Swift observation performed ~ 2 d before the 
beginning of the RXTE coverage. Assuming that the inner disc radius lies in between the neutron star radius and the corotation radius 
while the source shows pulsations, we estimate the magnetic field of the neutron star to be within ~ 2 x 10 s G and ~ 2.4 x 10 10 G. 
From this estimate, the value of the spin period and of the observed spin-up rate, we associate this source with the still poorly sampled 
population of slow, mildly recycled, accreting pulsars. 

Key words, stars: neutron — stars: rotation — pulsars: individual (IGR 117480-2446) — X-rays: binaries 



1. Introduction 

The dense environment of a globular cluster and th e resulting 
frequent stellar encounters (Mev lan & Heggiel [l997) make the 
production of binary systems hosting a compact object very effi- 
cient. Terzan 5 is proba bly one of the densest and metal-r ichest 
cluster in our Galaxy (ICohn et al]l2002tlOrtolani et al.ll2007l) . as 
clearly indicated by the large number of r otation-powered mil- 
lisecond pulsars dis covered there ( i 30, Ran som et al.l 12005: 
iHessels et al.l I2006I) . The cluster also contains at least 28 
discrete X-ray sources, a substantial number of which can 
be identified as qu iescent low-mass X-ray binaries (LMXB, 
iHeinke et al. 2006|). According to the recy cling scenario (see, 
e.g. iBhattacharva & van den Heuvei I1991I) . the population of 
rotation-powered millisecond pulsars and low-mass X-ray bi- 
naries (LMXB) share an evolutionary link, because the former 
are thought to be spun up by the accretion of mass and angular 
momentum in an LMXB. Accreting pulsars in LMXB are usu- 
ally found with periods clustering in two distinct groups. So far, 
13 sources have been discovered with spin periods lower than 
10 ms and were the refore termed accretion-powere d millisec- 
ond pulsars (see, e.g JWiinands & van der Ms|[l998). However, 
a smaller number of sources are found with longer periods and 
correspondingly higher estimates of the neutron star (NS) mag- 
netic field. 

So far, the only bright transient LMXB k nown in the clus- 
ter Te rzan 5 was the burster EXO 1745-248 (Makishim a et all 
1 1 98 lb . The first detection of a new outburst of a source in this 
cluster was made o n 2010 October 10.365 with INTEGRAL 
dBordas et al.l 120101) and was tentatively attributed to EXO 
1745-248. Follow-up Swift observations refined the source po- 
sition, and a comparison with the position of sources previ- 
ously known thanks to Chandra observations of the cluster 



suggested instead a different associatio n (IHeinke et al.1 120101: 
iKennea et ail 120101: iPoolev et all |2010|) . The X-ray transient 
is then considered a newly discovered source and named as 
IGR J17480-2446 (J17480 in the following). A coherent sig- 
nal at a period of 90.6 ms was detected thanks to observa- 
tions performed with the Rossi X-ray Timing Exp lorer (RXTE 
in the following, Strohmaver & Markwardt 2010). A signal at 
the same period i s present also during the several bursts that 
the source shows ( Altamiran o et al.1 1201 Obi) , while burst oscil- 
lations have never been observed from NS rotating at periods 
exceeding a few ms. A sudden decrease of the flux was tenta- 
tively attributed to an eclipse of the source by the companion 
(Strohmaver & Markwardt 2010). However, eclipses were not 
observed during subsequent observations, and the earlier flux de- 
creas e was identified with a lunar occultation (Strohmave r et al.l 
1201 Oi S10 in the following). 

Below we present the first analysis of the properties of the 
coherent signal emitted by this source, using RXTE and Swift ob- 
servations performed between 2010 October 10 and November 
6, and give a refined orbital and timi ng solution of the pulsar 
with respect to those first proposed (Papitt o et alj|2010l P10 in 
the following, S10). 



2. Observations and data analysis 

After the source discovery on 2010 October 10.365 
dBordas et al.l 1201 Ol) . RXTE started monitoring the source 
on MJD 55482.010 (October 13.010; aU the epochs reported in 
this paper are given with regard to the Barycentric Dynamical 
Time, TDB, system). We present an analysis of the observations 
performed until MJD 55506.359 (Obsld 95437-01-12-01), for 
a total exposure of 206 ks. In this time interval, a large number 
(> 300) of X-ray bursts are observed with a recurrence time 
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Fig. 1. J17480 in the time interval considered here as observed 
by the PCU2 of the PCA aboard RXTE. 

decreasing from > 26 to ~ 3 min as the X-ray flux increases. 
The analysis of the bursts shown by the source, as well as of 
its aperiodic timing properties, will appear in a companion 
paper. The 2.5-25 keV light curve recorded by PCU2 of the 
Proportional Counter Array (PCA) on-board RXTE, with the 
burst intervals removed and background subtracted, is plotted 
in FigJT] The count rate increases during the first days of the 
outburst, reaching a peak value of ~ 1700 s at MJD 55487.5, 
and then decreases to a value of ~ 800 s with an exponential 
decay time scale of ~ 5 d. Quite complex dipping-like structures 
also appear especially between MJD 55485 and 55490 as the 
flux shows sudden variations up to 75% on timescales of ~ 10 
min. 

The combined X-ray spectrum observed by the top layer of 
the PCU2 of the PCA (2.5-25 keV), and by Cluster A of the High 
Energy X-ray Timing Experiment (HEXTE, 22-50 ke V), can be 
well modelled by a sum of a blackbody an d a Comptonized com- 
ponen t, which we model with compps dPoutanen & Svenssonl 
119961) . Throughout the observations considered here, the spec- 
trum softens significantly as the source evolves towards higher 
luminosities. The unabsorbed total flux, extrapolated in the 0.1- 
100 keV band, rises from 0.47(3) x 10" 8 to a maximum level of 

I. 89(4) x 10" 8 erg cirT 2 s _1 , observed during the observation of 
MJD 55487.5. All the uncertainties on the fluxes given here are 
quoted at a 90% confidence level. 

To analyse the spin and orbital properties of the source we 
consider data taken by the PCA in Event ( 122jus temporal resolu- 
tion) and Good Xenon (lyus temporal resolution) packing modes. 
We discard 10s prior, and 100s after the onset of each type-I 
X-ray burst. The time series were also preliminarly barycentred 
with respect to the solar system barycentre using the available 
RXTE orbit files and assuming the best Chandra estimate of 
the source position, RA=17 A 48'" 4.831 s , DEC=-2 4° 46' 48.87" 
with a n error circle of 0. 06" (lcr confidence level, Heink eet alJ 
2006; iPoolev et alJbOlOl) . A coherent signal at a frequency of 

I I. 045(1) Hz [equivalent to a period of 90.539(8) ms] is clearly 
detected in the power spectrum at a Leahy-normalised power of 
=2 1.5 x 10 4 . A first orbital solution is obtained modelling the 
observed modulation caused by the orbital motion, 

f 2ttx 1 
P(t) = P < 1 + [cosra + ecos(2m -co)]> . (1) 

y Porb j 

Here Pq is the barycentric spin period of the source, x = a sin i/c 
the semi-major axis of the NS orbit, P or t, the orbital period, 
m = 2n(t - T*) I Port, the mean anomaly, T* the epoch at which 
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Fig. 2. Pulse profiles and best-fitting harmonic decomposition, 
obtained by folding observations performed by RXTE on MJD 
55482 (left scale, solid line) and MJD 55483 (right scale, dashed 
lines) around Pp = 90.539645 ms. Both profiles are evaluated 
after background subtraction and are normalised to the current 
average flux. 

m — 0, e the eccentricity and lo the longitude of the periastron 
measured from the ascending node. The periods P(f) are esti- 
mated by performing an epoch-folding search on 1.5 ks long 
data segments (for a total of 127) around the periodicity indi- 
cated by the po wer spectrum. The resulting variance profiles are 
fitted following lEeahvl (fl987h and the uncertainties affecting the 
period estimates are evaluated accordingly. The best-fitting or- 
bital solution we obtain with this technique is shown in the left- 
most column of Table Q] 

The time series were then corrected for the orbital motion 
with these parameters and folded around the best estimate of 
the spin period over 300s data segments (for a total of 717). 
The pulse profiles could generally be modelled using up to three 
harmonics. The pulsed fraction is observed to greatly vary in 
between the first (MJD 55482.010 to 55482.043) and the other 
RXTE observations. During the former the pulse fractional am- 
plitude of the first harmonic is very high [A\ 0.252(2)], while 
the second harmonic is detected at a much lower amplitude 
[A 2 = 0.016(2)]. In subsequent observations the amplitude of 
the first harmonic drastically decreases to values between 0.02 
and 0.04, whereas the second harmonic amplitude remains sta- 
ble and a third harmonic is sometimes requested by the profile 
modelling. To show this, we plot in Fig. [2] the pulse profiles 
calculated over the observations performed during MJD 55482 
(solid line) and MJD 55483 (dashed line), with the latter profile 
shown at a magnified scale. The pulsed fraction decrease is ev- 
ident, as are the variation of the shape of the peak at rotational 
phase ~ 0.75. 

To increase the accuracy of our timing solution we model 
the temporal evolution of the phases evaluated on the first and 
second harmonic of the pulse profiles with the relation: 

<f>(t) = + (v - v F ) (t ~ Tref) + -< v >(t - T ref f + R or b(t).(2) 

Here T re f is the reference epoch for the timing solution, vo is the 
pulsar frequency at the reference epoch, Vf = I/Pf is the fold- 
ing frequency and < v > is the average spin frequency derivative. 
The term R or b{t) describes the phase residuals induced by a dif- 
ference between the actual orbital parameters, namely x, P or b, 
T*, e sin a> and e cos co, and those used to correct the time series 
(see e.g. lDeeter et al.lfl98ll) . Once a new set of orbital parame- 
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Table 1. Spin and orbital parameters of IGR J17480-2446. 
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<v> (lQ- i2 Hz/s) 


< 16 


1.22(1) 


1.68(1) 


< 




asin i/c (lt-s) 


2.498(5) 


2.4967(3) 


2.4973(2) 




-3 


Porb (hr) 


21.2744(8) 


21.2745(1) 


21.27454(8) 






T*- 55481.0 (MJD) 


0.7805(4) 


0.78033(6) 


0.78048(4) 







e 


< 0.02 


< 7 x 10~ 4 


< 6 x 10~ 4 






f(M„M 2 ,i) (M ) 


0.0213(2) 


0.0212587(8) 


0.021275(5) 




-1 


x l /dof 


156/121 


7879/662 


3287/566 


o 
<J 


-2 



Notes. Numbers in parentheses are lcr errors on the last significant 
digit. Upper limits are evaluated at 3<x confidence level. The uncertain- 
ties of the timing parameters have been scaled by a factor -\/xJ to take 
into account a reduced^ 2 of the best-fitting model larger than 1. The or- 
bital solution obtained modelling the Doppler shifts on the pulse period 
is given in the leftmost column, while the best timing solution evalu- 
ated from the phase evolution of the first and second harmonic are given 
in the central and in the rightmost column, respectively. The reference 
epoch for these timing solutions is MJD 55483.186. 
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Fig. 3. Evolution of the first (top) and second (bottom) har- 
monic phases obtained by folding the orbitally corrected time 
series around Pf = 90.539645 ms (see text). The solid lines 
show the best-fit models from Eq. (0. 



ters is found, it is used to correct the time series, and the resulting 
phases are modelled again with Eq.©. This procedure is iterated 
until the phase residuals are normally distributed around zero. 

Although no residual modulation at the orbital period is ob- 
served, the phases of the first harmonic are strongly affected by 
timing noise. The reduced x 2 we obtain modelling their evo- 
lution with Eq.© is extremely large (^ 11.9 over 662 d.o.f.). 
Such a behaviour is most probably caused by to pulse shape 
changes like the one shown in Fig [2] The second harmonic 
phases appear to be less affected by timing noise, resulting in a 
reduced^ = 5.8 (566 d.o.f.). We argue that the second harmonic 
phases are better fitted with respect to the first harmonic because 
of the greater stability of this component (as alre ady observed 
in some accreting millis econd pulsars, see, e.g., Burderi et al.l 
120061 [Riggio et al. 2008). The best-fitting parameters calculated 
over the phase evolution of the first and second harmonic are 
quoted in the central and rightmost column of Table [T] In Fig. 
[3] we show the phases of both harmonics, when the observa- 
tions corrected for the orbital motion of the source are folded 
around Pf = 90.539645 ms. The phase evolution is clearly 
driven by at least a quadratic component. A consequence of tim- 
ing noise is that the spin frequency and its derivative, estimated 
over the two harmonic components, are significantly different. 
We quote conservatively a value of vo = 11.0448854(2) Hz 
that overlaps both frequency estimates, and use a spin frequency 
derivative between 1.2-1. 7xl0~ 12 Hz s _1 in the discussion be- 
low. However, it is worthwhile to note that the orbital parame- 
ters are entirely consistent between the two harmonic solutions, 
which supports the reliability of these estimates. The solution we 
obtained is entirely compatible with, but more precise than, those 
proposed by PlQjand S10. Given the accuracy of the source po- 
sition considered here (0.06"), the systematic uncertainties in- 
troduced by the position indetermination on the measured values 
of spin frequency and of its derivative (e.g. lBurderi et al.|l2007l) 
are cr pos v i 3 x 10~ 10 Hz and cr pos ;j , i6x 10~ 17 Hz s _1 , respec- 
tively. Finally, as the cluster moves towards the solar system at 
a velocity of 85 + 10 km s _I ([Ferraro et alJ l2009). the measured 
value of the spin frequency is affected by a systematic offset of 



1 There is an offset between the values of frequency and epoch of 
mean longitude quoted by P10 and those presented here, as theirs were 
not referred to the TDB reference system. 



— 1-3 x 10~ 3 Hz, though this is unimportant when making con- 
clusions about the source properties. 

In order to extend the range of fluxes at which the source 
was observed and pulsations were detected, we also analysed 
three Swift observations (Obs. 00031841002, 00031841003 and 
0003 1 84 1 004) in which the XRT observed in Windowed Timing 
(WT) mode, with a temporal resolution of 1 .7ms. The Swift XRT 
started monitoring the source on MJD 55479.737, more than 
two days before RXTE, with the first observation in WT mode 
starting at MJD 55479.802 for 2 ks. After applying barycentric 
corrections for the satellite orbit, then correcting for the source 
orbital motion and selecting photons from a 50 pixel wide box 
around the source position, a pulsation is clearly detected at a 
period of P$ = 90.5395(2) ms by means of an epoch-folding 
search. The XRT signal is particularly strong and consistent with 
that seen by RXTE during its first observation, with a pulse pro- 
file modelled by a sinusoid of amplitude 0.23(2). Pulsations were 
also searched for in the subsequent ~ 1 ks long XRT WT ob- 
servations starting on MJD 55484.767 and 55485.357. Only a 
weak signal was detected in the latter at a fractional amplitude 
of 0.018(4), still compatible with that seen by RXTE at those 
later times. 



3. Discussion and conclusions 

We reported on the spin and orbital properties of the newly dis- 
covered accreting pulsar, IGR J17480-2446. Its 90.6 ms period 
makes it the first confirmed accreting pulsar in the range 10- 
100 ms. Pulsations were detected in all observations performed 
by RXTE, as well as in two out of the three Swift observations 
performed in WT mode presented here. The pulsed fraction is 
observed to drastically change on a timescale of i 1 d, after 
the RXTE observation performed on MJD 55482. While previ- 
ously both Swift and RXTE observations revealed a strong signal 
dominated by a first harmonic component of fractional ampli- 
tude as large as 0.25, later observations at higher fluxes per- 
formed by both satellites never detected an amplitude £ 0.03. 
Simultaneously the pulse shape changes and becomes more 
complex. This behaviour is suggestive of a change of the geo- 
metrical properties of the flow in the accretion columns above 
the NS hot spots. 
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Because pulsations are detected throughout the observations 
shown here, an estimate of the NS magnetic field strength can be 
made. For accretion to proceed and pulsations to appear in the 
X-ray light curve, the inner disc radius, /?,„, has to lie in between 
the NS radius, Rns, and the corotation radius, Rc, defined as the 
distance from the NS at which the velocity of the magnetosphere 
equals the Keplerian velocity of the matter in the disc. For a 
larger accretion radius, accretion would be inhibited or severely 
reduced by the onset of a centrifugal barrier. For a pulsar spin- 
ning at 90.6 ms, R c = (GMP 2 /4tt 2 ) 1/3 = 33&m\ / l km, where 
mi a is the NS mass in units of 1.4 M Q . Defining the inner disc 
radius in terms of the pressure equilibrium between the disc and 

1/7 —2/7 —2/7 4/7 

the magnetosphere, one obtains =* 160 m 14 R 6 L 37 // 2 g 
km dBurderi et alJ uOOl ). where R(, is the NS radius in units of 
10 km, L37 the accretion luminosity in units of 10 37 erg s _1 , 
and yU28 the magnetic dipole moment of the NS in units of 10 28 
G cm 3 . Extrapolating the fluxes observed by RXTE to the 0.1- 
100 keV ba nd, and assuming a distance of d = 5.5 ± 0.9 kpc 
to Terzan 5 dOrtolani et al]l2007D . we estimate a maximum and 
minimum bolometric luminosities of 1.7(1) x 10 37 d? 5 erg s" 1 
and 6.8(1) x 10 37 d 2 erg s _1 , during the time covered by the 
observations considered here. Assuming that the X-ray lumi- 
nosity is a good tracer of the accretion power and imposing 
Rns < Rin ~ Rc, we obtain 



0.02 m~2 H R 6 



1/4*9/4. 



4.8 m\^Rl /2 d 5 



(3) 



The upper limit on the magnetic dipole can be reduced con- 
sidering that pulsations are detected also in a Swift observa- 
tion taking place ~ 2 d earlier than the first RXTE observation. 
iBozzo et al.1 (12010) estimated the source flux in that observation 
as 4.5(2) x 10~ 10 erg cirT 2 s~' (1-10 keV). This value is a factor 
~ 4 lower than the value obtained extrapolating the spectrum of 
the first RXTE observation in the same energy band. Assuming 
that this ratio holds also for the bolometric luminosity of the 
source, we get to an upper limit on the magnetic dipole moment 
of ^ 2.4 x 10 28 G cm" 3 . The limits thus derived translate to a 
magnetic surface flux density between ~ 2x 10 8 and ~ 2.4 x 10 10 
G. The upper bound of this interval can be overestimated be- 
cause the exact flux at which the pulsations appeared is unknown 
at present. Monitoring the presence of coherent pulsations as a 
function of the flux when the s ource fades will probably allow 
us to derive a tighter constraint. Alta mirano et a D(H)10a) have 
also reported the presence of a kHz QPO at ~ 815 Hz (10-50 
keV) during the RXTE observations performed on MJD 55487. 
Under the hypothesis that this feature originates in the inner- 
most part of the accretion disc, it indicates an inner disc radius 
Ri„ i 20 m^ 3 km. As the luminosity we estimated during that 
day is 6.8(l)x 10 37 d 2 5 erg s _1 , this would imply a magnetic field 
i7x 10 8 ds 5 G if the disc is truncated at the magnetospheric 
radius. 

Despite the presence of timing noise, the analysis of the 
phase evolution over the ~ 24 d time interval presented here 
clearly indicates the need for a quadratic component to model 
these phases. Interpreting this component as a tracer of the NS 
spin evolution, we thus conclude that the source spins up while 
accreting. Values of the spin-up rate between 1.2 and 1.7xl0~ 12 
Hz s" 1 are found, depending on the harmonic considered. This 
discrepancy is probably due to the effect of timing noise. These 
values are compatible with those expected for a NS accreting 
the Keplerian disc matter angular momentum given the observed 
luminosity, v ^ 1.5 x 10~ 12 (L 37 /5) (7?,„/70/fcm) 1/2 R 6 m L4 
Hz s _1 . Here 745 is the NS moment of inertia in units of 10 45 



g cm . The observed spin period and the magnetic field we 
estimated place this source between the population of "classi- 
cal" (B > 10 n G, P > 0.1 s) and millisecond (B a 10 8 -10 9 G, 
P - 1.5-10ms) rotation-powered pulsars. The observation of 
a significant spin-up at rates compatible with those predicted 
by the recycling scenario further supports the identification of 
this source as a slow, mildly recycled pulsar. We note that the 
only other two accreting pulsars with similar, though signifi- 
cantly different param eters, are GRO J 1744-28 (Ps = 467ms, 
B ^ 2.4 x 1 11 G.ICuill 19971). and 2A 1822-371 (P s = 590 ms, 
B <10 U G. Uonker & van der Klisll200l) . 

The orbital parameters we measured for the NS in J 17480 al- 
low us to derive constraints on the nature of its companion star. 
With a mass function of f(Mz',Mi,i) - 0.02 M , a minimum 
mass for the companion can be estimated to be as low as 0.41 
M Q for an inclination of 90°, and an NS mass of 1.4 M . Since 
the source shows no eclipses, the inclination is most probably 
^ 80°, and the lower limit increases to 1112 = 0.16 + 0.26mi 4, 
where ni2 is the mass of the companion star in solar units. An up- 
per limit can be obtained if the companion star is ass umed not 
to ove rfill its Roche lobe. Using the relation given by Eggleton 
(119831) and the third Kepler law to relate the Roche Lobe ra- 
dius to the orbital period and to the companion mass, R12 - 
O.55(GM ) 1/3 (/ J ori /27r) 2/3 m 1 .4 ^ 2/3 (l + q) 1/3 /[Q.6q 2 ' 3 + log(l + 
q 1 ^ 3 )], where q = M2/M1, and assuming the companion star fol- 
lows a main sequence mass-radius relation, R2/R& ~ (M2/M0), 
yields a maximum mass of 2.75 M for the companion when 
mi 4 = 1. This upper limit is indeed higher than the maximum 
mass expected for a main sequen ce star belonging to one of the 
two stellar populations found by Ferra ro et al.l d2009h in Terzan 
5. One has in fact 1112 ~ 0.95 if the companion of J 1 7480 belongs 
to the older population (t=10 Gvr. ID' Antona et al.ll2010l) . while 
m.% i 1.2 and ntz i 1 .5 if it belongs to a younger population of 6 
and 4 Gyr, respectively (D'Antona, priv. comm.). We conclude 
that a reasonable upper limit for the companion-star mass is 1.5 
M Q , possibly a main sequence or a slightly evolved star. 
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